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TECHNICAL MEMORANDUM 


ANALYSIS OF NITROGEN CONDENSATION 
IN AN EXPANDING NOZZLE FLOW 

INTRODUCTION 


In a recent repoi'l [1], experimental results of vapor-liquid condensation 
in small, contoured no/zles were presented. In these exiieriments, condensation 
takes place very close to the saturation point, leaving almost no supersatui’ation 
in the How. Thi-s is in contrast to most of the experimental data reported else- 
where |2-liJ]. To assess the factors which caused the discrepancies, a How 
analysis based on the condensation model of Griffin and Sherman [ 11] has 
been pcrfoianed. The discrepancies between analysis and experiment indicate 
that there are factors which are not c-onsidered in the analytical model. 

The two factors known to have a major effect on the flow condensation 
are the rate of expansion and the decree of impurities. The former is a fluid 
d\mamic factor and can be accounted for analytically usin^ the equations of fluid 
flow. Usinn experimental data obtained primarily from supersonic and hyper- 
sonic wind tunnels, Daum and Gyarmathy Ilf)] oljtained a correlation between 
the onset of condensation and the rate of expansion of the nozzle. The exi)ansion 
rate also affects the dro])let p,rowth onc-e the onset of condensation has been 
reached. The second factor is a thermodynamic phcnomenoti and c an l)c explained 
by the kinetic theory of liqidds lldj. Experimental studies conducted l)V Arthur 
and Na^amatsu joj indicate that the foreign particles act like additional nuclei 
which enhance condensation. Eor llcnvs with lar^e amounts of foreign particles, 
condensation may become so rapid that the How follows the vapor pi'cssure 
curve, leaving no supersatuiation at all. 

The present condensation model usc-s the thcoi v of micleation with the 
steady one-dimensional equations of motiem for diabatic Hows. Computer 
analyses for Hows in oiu'-dimensional nozzles usini; these c'quations can he 
found in other re])orts j 17-Ihj . 



ANALYSIS 


The analysis uses the equations of one-climcnsif)nal diabatic flow, the 
theory of nueleation, and the droplet nn)wth rate equation. The basic assump- 
tions used are as follows: ( 1) the flow is steady, one-dimensional, and the 
noz/.le wall is frictionless and non-heat eonductinK; (2) the vaiK)r is a perfect 
gas with known thermodynamic constants; (2) the flow is isentropic before the 
onset of condensation, and diabatic thereafter; (4) the condensed mass is li(|uid, 
and its interaction with the flow is negligible; ( li) tlie volume of the condensed 
jjhase is negligible comi)ared with the U)tal volume. With these assumptions, 
the How equations can be written as 1 17|: 
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where p, T, 1', A, and p :ire flow quantities measured along the streamtube 
under consideration; L is the latent heat of evaporation; and g is the mass 
fraction in the condensed phase. 

i:(|u:itions (l) through (a) serve as the isentro|)ie condition when g and 
dg are zi'i'o, :md the dialjatic ef>ndition when the condensed phase is i)resent 
(g ^ 0). Iteplacing tlie differenti:ils in these equations by differinec's and 
rearranging, the following equations for the step ehange of |)ressure :ind tempi-ra- 
ture can bi- obtained: 



and 



(•i) 




(■•*) 


are the fractional de\iation of pressure and temperature frnm tluir coi n spond- 
isentnipic conditions, rrom the inc-asured flow guantitii s alonu the noz/lc , 
equations ((!) through ( '') ^i'e tlie valuc-s of . and . From equations ('') 

and (‘J) it is obvious that ‘ ^ ' I f*’*' ^ ^ T M 1 • llc nct , 

in a suiiersonic flow the- temperature* measurements will Kive better indications 
of onset of condc*nsation. 

To calculate the inc*rement AUi the equations for the* rate of nucleatir)n 
and the droplet ^;rowtli rate* aie used. Fiiom the* Kine tic the*oiy descniilion of 
Volmer ;> condensation nucleus is formed s|K)ntaneously upon collisions 

lH'twee*n clusters of vapor mole*cule*s. U lu*n tlie* nue*leus, so formed, re aches 
a critical radius, it will continue to ki’ow upon fiirtlier collisions with \apor 
malecule-s; othe*neise, it will e*vai)orate and diminisli eve ntuall\ . .Assuming that 
( 1) the collision pndialiility on a unit surface is ^ii\e*n l)\ an eepiilibi'ium pailie h 
elistril)utie»n, (2) tlie saturation vapor pressure ee|uals that of a drople I ol inlinili 
radius, and (.‘5) the radius of the di^tplct is smalle r tlian the* me*an fn e* jiath ol 
the* vapor, the following; eq iati<)ns can be* obtaine*el [ ]7|: 
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where Jj is the number of critical drops formed per unit volume per unit time, 

k the Boltzmann constant, N the Avadro's number, p the molecular wcinht of 

J\ 

the vapor, a the surface tension of the liquid, and r^ is the critical drop 
radius. 

The number of droplets formed in a step increment of Ax, is then ^i'en 
by 


N. = J.A.Ax, 
1 1 i i 


( 12 ) 


By a consideration of heat transfer in a rarefied fjas, the ^^^o^^lh rate of a drop- 
let can be found as 



( 13 ) 


where O' is an accommodation factor and T is the temperature of the liquid. 

s 

By summing up the j;rowth in the step increment Ax^ of droidcts formed 

in all the previous steps, the net mass fraction condensed in this step increment 
can be obtained as 
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The surface tension, a , used in c<iuations (10) and (11) is a fimetion of 
the droplet radius. A relation niven by Tolman 120] is usetl here; 
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where a is the surface tension when the drop radius api)roachcs infinity. The 

GO 

parameter D is the Tolman constant which is closely related to the intermolecu- 
lar distance in the liciuid dmp. The value of the order of 10“'* cm has Ix'cn 
suRfjested in Hcfcrcnces 21 and 22. However, the calculations used arc not 
sufficiently reliable b) give an exact value. Variations in the value of Tol man’s 
constant are thus included in this analysis. 


To start the calculation, the rate of expansion (distribution of area A) 
of an cKjuivalent no/./.le is obtained from the Mach number distribution of the 
actual nozzle flow. Using this area distribution as an input, and with a given 
chamber condition, isenti’opic relations can In- used until saturation condition 
is reached. The nuclcation relations, equations (10) and (11), are then used 
to calculate the initial condensation at each increment f)f x thereafter. Hy com- 
I)aring the dn)plet size with a given critical drop size, the |K)int <»f onset of cf)n- 
densation can also be defined. The dn)|)let growth equations (12) and (10) are 
then included in the calculation |)ro\1ding complete flow |)arameters at each 
step Ax. 


The effect of foreign nuclei on the onset of condensation can be included 
by assuming that all the nuclei arc of the same size, and an equivalent radius 
r^j for the nuclei can be defined. With a given density and mass fracti<<n 

7 | for the foreign nuclei in a flow t)f mass flux m, this equivalent radius can be 

related to the number flux N , and surface flux S , of the nuclei by the follf>wing 

d d 

etpiations: 
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Note that the factor lui; 


gives the volume Ilux of Liu iori'ign luielei. 


liy ni'ulcctinn the Interaction between the foreign nuclei and the vajKjr 
flow, the only effect caused by the presence of the foreign nuclei is Ui pnivide 
additional surfaces on which the vajxjr can condense. I’nder this circumstance, 
the vapor starts to condense on the forelj^i nuclei Immediately uix>n reachinji 
the saturation condition. Homoneneous condensation may still occur later ard 
the same proceduix' mentioned in the j)ie\1ous paragraph will follow. 


RESULTS AND DISCUSSIONS 


Analysis of nitrf)^^en flow in the contoured no/./Ie used in Hi ferenee 1 has 
been |)erlbrnu’d for two sets of chamber conditions: ( 1) - (isi. 'j N/em‘ 

('jsp psia), lisu K (noPIt), and (li) I»^ = Gls. l N/cm^ j»sia), ^ 

27s. K (nuPR). These conditions correspond to 'I'est Puns 9 and P» in Itefer- 
ence 1 (Table 1). Pesults of j)arametric studies for v.arious values of the mass 
accommodation factor, n , and the Tolman's e*»nstant, D, are also presenti-d. 


The distriljution of temperature, jiercentaue condensed, and tlu- rate of 
nucleation for clianil)er condition ( l) are j’iu-n in Piuiires I, 2, and ii. It is 
evident fn)!!! these fi'iures that the onset of condensation moves upstream as the 
value of Tolman's ctuistant increases. An upper Imuod of approximately I) 

•l.n ■ ]<)~^ cm, however, ck>es exist, beyond which the nunu'rical scheme lireaKs 
doATJ. This m:iy be due to tlie laiTe amount of condensati(>n which occurs in a 
single stepsi/.e that stops the calculation. 


The effect on thv onset of condensation due to changes in tlu* aceonunoda- 
tion bictor is rather small. Its main effect, rather, shows in the laiti- of iuule:i- 
tion. I’n)m I'inure 2 it can l)c shown Unit for laiTc accommod.ition fachirs tlu- 
nucleation occurs at hinh rates over :i veiy luirnjw re;iion, :> suiidi-n and 

drastic change in the thermodvmimic states of tlie lluid. At very low \ tilues, 
e.^;. , n =<•.()!, tlu rate of nucleation beconu‘s so snuill that cS i ffi’ct on tin 
flow is com|)ensated by the rate of expansion in the How. ('ontiiuious nuc-k*ati(»n 
at moderate ••ates cxtenils all ihe way toward tlr* i-xit of the no/./le. I'his 
ohvnomenon can :ilso be shown clearly in the pivssnri--tempei:ituri- (list libulion 
^,iven in Ti^;ure 1 where :dl the curves a|)proach a parallel to the s:iturati(m 
curve. The tnmsition from the is( titn>pe to the parallel is where the thenno- 
dynamic n(*n-e(juilibriiim occui’s. The diminishing; of this tr:insition :e;;ion loi' 
small (. 's indic:ites a state of thei niodynamic (|u;isi-e(|idl iliriuin with super- 
saturation Kejit :it :i const:int level. 


(1 



The results presented in Figures 1 through 4 for the chamber condition 
( 1) are typical of all chamber conditions. Similar results have been obtained 
for chamber condition (2), but are not presented in this paper. 

Analysis of heterogeneous condensation with foreign nuclei are per- 
formed under six different conditions. These conditions, given in Table 2, cover 
a wide range of number density of particles and of particle size distribution. 

The resulting static pressure distribution for chamber condition { 2) with a = 

0. 75 and D = 4. 0 x 10“^ is given in Figure 7. It is obvious that under such a 
chamber condition, the condensation is affected very little by the variations in 
the foreign particles. 

Experimental measurements of the static pressure distribution from 
Reference 1 are presented in Figures 5, G, and 7 and compared with analytical 
results. Results indicate that the flow condensation occurs sooner than predicted 
by the analysis. This difference cannot be accounted for even if foreign particles 
were introduced in the analysis (Fig. 7). It is thus concluded that certain 
phenomena other than the spontaneous condensation of nitrogen exist in these 
flows. 


The tests conducted in Reference 1 have chamber conditions which lead 
to high saturation temperatures and pressures. A quick calculation for an 
adiabatic expansion from chamber condition (l) yields a saturation Mach number 
of 3. 8 and Reynolds number per unit length of 3 x 10“® cm“^. For a nozzle with 
a length scale of 1 cm as the one used in Reference 1, such a combination 
requires the flow to be treated as a continuum [23). This indicates that the 
free molecular assumptions used in the analytical model are no longer valid 
under this circumstance. The addition of viscous effect and molecular inter- 
action arc thus needed to better describe the flow. This observation is by no 
means an isolated one, as similar deviation from analysis (Fig. 8) is also 
found in the data of Goglia and Van W'ylcn [7J. Altliough the deviation was 
explained in Reference 15 as due to the possible presence of foreign particles, 
it could well l)c the ivsult of ncglcciing Llie effect of \1scosity imd molecular 
interaction. 



TABLE 1. EXPERIMENTAL CONDITIONS HI 


Run No. 

Vapor 

P 

c 

N/cm^ (psia) 

T 

c 

K (° R) 

9/1(1) 

Nz 

G81.9 (989) 

280.0 (504) 

10/0(8) 

Nz 

548. 1 (795) 

278.3 (501) 

14/0(10) 

Nz 

475.7 (G90) 

281.7 (507) 

11/0(8) 

Nz 

413.7 (GOO) 

280.0 (504) 

13/0(10) 

Nz 

344.7 (500) 

280.0 (504) 

12/0(5) 

Nz 

275.8 (400) 

279.4 (503) 


TABLE 2. FOREIGN NUCLEI CONDITIONS 


Case 

7? ( " ) 

(cm) 

drcps 

N . 

d s 

2 

• cm^ 
S . 
d s 

1 

0 

— 

— 

— 

2 

0.01 

5 X 10"'^ 

9. 2 X 10*^ 

2. 89 X 10® 

3 

0.01 

5x ]0“” 

9. 2x 10* 

2. 89 X lo'^ 

4 

1 

5 X 10“'’’ 

9. 2 X 10‘® 

2. 89 X lo'^ 

5 

1 

5x 10"^ 

9. 2x 10*9 

2. 89 X 10® 

() 

1 

5x 10"® 

9. 2x K)22 



2.89X 10® 
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I'iuure 7. Variation of static pressure distribution with foreijjn particles (run lO) 
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